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Glueballs, strings and topology in SU(N) gauge theory 

M. Teper*" * 

'^Theoretical Physics, University of Oxford, 

1 Keble Road, Oxford 0X1 3NP, United Kingdom 

I show how one can use lattice methods to calculate various continuum properties of SU(A'^) gauge theories; 
in part to explore old ideas that N = 3 might be close to X = oo. I describe calculations of the low-lying 
‘gluebair mass spectrum, of the string tensions of fc-strings and of topological fluctuations for 2 < N < 5. We 
find that mass ratios appear to show a rapid approach to the large-A^ limit, and, indeed, can be described all 
the way down to SU(2) using just a leading 0{1/N^) correction. We confirm that the smooth large-X limit we 
find is confining and is obtained by keeping a constant ’t Hooft coupling. We find that the ratio of the k = 2 
string tension to the k = 1 fundamental string tension is much less than the naive (unbound) value of 2 and is 
considerably greater than the naive bag model prediction; in fact it is consistent, within quite small errors, with 
either the M(-theory)QCD-inspired conjecture that au oc s\n{-Kk/N) or with ‘Casimir scaling’. Finally I describe 
calculations of the topological charge of the gauge fields. We observe that, as expected, the density of small-size 
instantons vanishes rapidly as N increases, while the topological susceptibility appears to have a non-zero N = oo 
limit. 


1. INTRODUCTION 

How SU(N) gauge theories approach their N = 
oo limit and what that limit is, are interest¬ 
ing questions ||^, in themselves, whose answers 
would, in addition, represent a significant step to¬ 
wards addressing the same question in the context 
of QCD. Accurate lattice calculations in 2-1-1 di¬ 
mensions Q show that in that case the approach 
is remarkably precocious in that even N = 2 \s 
close to N = oo. Such calculations have to be 
very accurate because for each value of N one 
has to perform a continuum extrapolation of var¬ 
ious mass ratios and then these are compared and 
extrapolated to iV = oo. Earlier D=3-|-l calcula¬ 
tions 1^,^ were much too rough for this purpose 
even if their message was optimistic. Recently, 
however, this situation has been rapidly improv¬ 
ing and here I want to describe some of the things 
that one has already learned. 

Analyses of Feynman diagrams to all orders im¬ 
ply [Q, that the SU(N) theory will have a smooth 
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limit if we keep fixed the’t Hooft coupling 

A = g^N. (1) 

We would like a non-perturbative confirmation 
of this expectation. Again, we would like a 
confirmation of the assumption that the large- 
N theory remains confining, since much of the 
phenomenological argument that QCD might be 
close to V = oo relies on this being so Q. Since 
quarks are in the fundamental representation, the 
leading correction to V = oo is expected to be 
0{1/N). In the pure gauge theory it is expected 
to be 0(1/V^). (Again an expectation that one 
would wish to test non-perturbatively.) Thus if 
QCD is close to V = oo one would expect the 
SU(3) gauge theory to be close to SU(oo) as well. 
We can test this by calculating dimensionless ra¬ 
tios of physical masses as a function of N. There 
are also arguments and speculations (e.g. [||) 

that topological fluctuations are very different at 
large N and at small N ; this too we would like to 
test in lattice simulations. 

The calculations ||Q on glueballs, topology 
and fc-strings which I will describe have been per¬ 
formed in collaboration with Biagio Lucini. The 
strategy is simple: we calculate the relevant (con- 
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tinuum) properties of SU(2), SU(3), SU(4) and 
SU(5) theories and directly compare them. Since 
= 5 is large in the sense that the leading cor¬ 
rection oc 1/A^^ = 1/25 <C 1, it is not surprising 
that such a calculation turns out to suffice for 
our purposes. These calculations were intended 
as exploratory. A much larger calculation is now 
underway and this will provide much more infor¬ 
mation on the mass spectrum and its dependence 
on the number of colours. 


2. FROM LATTICE TO PHYSICS 


Our lattice calculations are entirely standard 
1^,^. The Euclidean space-time lattice is hyper- 
cubic with periodic boundary conditions.The de¬ 
grees of freedom are SU(iV) matrices residing on 
the links of the lattice. In the partition function 
the fields are weighted with exp{S'} where S is 
the standard plaquette action 

^ = (2) 

p \ / 


in which Up is the ordered product of the matrices 
on the boundary of the plaquette p. For smooth 
fields this action reduces to the usual continuum 
action with j3 = 2N/. For rough fields on a 
lattice of spacing a we can define a running lattice 
coupling gL{o) which reduces in the continuum 
limit to a coupling g[a) in our favourite scheme: 


2N a^o 2A 
^ - glia) ^ g^{a) 


( 3 ) 


Path Integrals using standard Monte Carlo tech¬ 
niques. On the lattice t = na and so we obtain 
the energies from eqn(^ as aE^ i.e. in units of 
the lattice spacing. In practice one needs to use 
several carefully chosen operators, and a varia¬ 
tional calculation . Since our volume is finite 
we must make sure that finite volume corrections 
are negligible. For the calculations I describe here 
such checks have been made. 

Having calculated some masses arui at a fixed 
value of a (i.e. at a fixed value of /3) we can re¬ 
move lattice units by taking ratios: arrii/amj = 
mi/rrij. This ratio differs from the desired contin¬ 
uum value by lattice corrections. For our action 
the functional form of the leading correction is 
known to be O(a^). Thus for small enough a we 
can extrapolate our calculated mass values 


rriiia) mj(0) 


)(a) 


,( 0 ) 


ca^mlia) 


( 5 ) 


where using different will make differences at 
O(a^) (as does the a-dependence of mfe(a)). We 
shall later show some explicit examples of such 
extrapolations. 

Having calculated such mass ratios in various 
continuum SU(fV) gauge theories, we can attempt 
to extrapolate to TV = oo using the fact that the 
leading correction is expected to be 0(l/iV^): 


TO, 


N 


7V2- 


( 6 ) 


If such fits are good for N > Nq then we may 
regard SU(Ao) as being close to SU(oo). 


Thus by varying the inverse lattice coupling f3 we 
can vary the lattice spacing a. 

To calculate a mass we construct some operator 
(pit) with the quantum numbers of the state (typi¬ 
cally this will be a linear combination of the prod¬ 
ucts of link matrices around some closed loops) 
and then use the standard decomposition of the 
Euclidean correlator in terms of energy eigen¬ 
states 

dt) = {p\t)piU)) = ^ |(H|^|n)pexp{-£;„t}(4) 

n 

where \n) are the energy eigenstates, with En the 
corresponding energies, and |H) is the vacuum 
state. We evaluate the corresponding Feynman 


3. GLUEBALLS 

In this exploratory calculation we focus on 
what are expected (from previous calculations in 
SU(2) and SU(3) |§,0) to be the lightest states: 
the lightest = 0++ and 2++ glueballs. We 
also calculate the mass of the first excited scalar 
glueball, which we shall refer to as the O'*'’*'*. In 
addition to these glueball masses we calculate the 
string tension, o^cr, of the flux tube between static 
sources in the fundamental representation (see 
the next Section.). This last is our most accu¬ 
rately calculated quantity and so we use it to form 
the dimensionless mass ratios ml ^a which we ex¬ 
trapolate to the continuum limit using eqn(^. In 
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Figure 1. The mass of the lightest scalar glueball, 
mo++, expressed in units of the string tension, cr, 
is plotted against the latter in lattice units. The 
a ^ 0 continuum extrapolation, using a leading 
lattice correction, is shown. 


Figure 2. Continuum scalar, tensor and excited 
scalar masses expressed in units of the string ten¬ 
sion and plotted against l/N'^. Linear extrapola¬ 
tions to TV = oo are shown in each case. 


Fig.|| I show you how this works for the lightest 
scalar glueball for N = 2,3,4, 5. 

While the mass of the scalar glueball is the 
most accurately calculated, because it is the light¬ 
est, it is also the one with (by far) the largest lat¬ 
tice corrections. Despite this it is already clear 
from Fig.|^ that for small a there is very little N- 
dependence for N > 3. To be quantitative we 
perform continuum extrapolations using eqn(p|), 
as shown in Fig.^ The results are shown in Fig]^ 
where we also show the extrapolations to N = oo 
using eqn(|). 

We see from Fig.|| that, as far as the lightest 
glueballs are concerned, all SU(iV) theories can 
be described by a modest leading 0{1/N‘^) cor¬ 
rection to the SU(oo) limit. 

We also see from these plots that the ratio of 
the glueball masses to the (square root of) the 
string tension has a finite non-zero limit as N 
oo. This tells us that the confining string tension 
remains finite and non-zero as N —s- oo. (Caution: 
numerical calculations only test confinement to 
some finite distance - which in the case of our 
higher-calculations is not yet very large.) 


4. ’T HOOFT COUPLING 


The analysis of diagrams [|p suggests that the 
way to achieve a smooth limit as TV ^ oo is by 
keeping the’t Hooft coupling, A = g^N, constant. 
Since the coupling runs, we should say that what 
we keep fixed is the running’t Hooft coupling, as 
defined on some scale I that is fixed in units of 
some quantity that partakes of the smooth large- 
N limit, such as the string tension. To do this we 
use eqn(||) which tells is that a suitable defintion 
of a running’t Hooft coupling is 


^/(a) = g]{a)N = 


2N^ 

0{ReTrUp/N) 


( 7 ) 


The extra factor involving the plaquette is a 
mean-field (or tadpole) improved version of (3 and 
the naive A(a) we would derive from it [pd| . It is 
necessary ||^ because the naive lattice coupling 
is known to be very poor in the sense of having 
very large higher order corrections. 

We extract A/ (a) for each of our various lattice 
calculations at various N and a . Our diagram¬ 
matic expectation is that if we plot a^a against 
A/(a) then, for large enough N, the calculated 
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Figure 3. The square root of the string tension 
in lattice units, a^/a, plotted against the T Hooft 
coupling, A/ = gjN. 


values should fall on a universal curve. In Fig.| 
we plot all our calculated values of a^ya, against 
the value of A/(a). We observe that, within cor¬ 
rections which are small except at the largest 
value of a, we do indeed see a universal curve, and 
so the diagrammatic expectation is supported by 
our non-perturbative calculations Q. 

5. STRINGS 

In a confining SU(N) gauge theory, the poten¬ 
tial between two static charges in the fundamen¬ 
tal representation is expected to grow linearly at 
large separations as V(r) = ar where a is the 
fundamental string tension. If we have charges in 
other representations then, in SU(3), they can be 
screened to either the fundamental or the trivial 
representation by gluons from the vacuum. So the 
fundamental string is the only stable one (assum¬ 
ing it to have the lowest tension). In S\J{N > 4) 
there are other stable strings. It is convenient to 
label a charge by the integer k if its wave func¬ 
tion transforms by a factor of under a global 
gauge transformation z belonging to the centre 
Zn of the group. Because gluons transform triv¬ 
ially under the centre, they cannot screen a charge 


of k into a charge of k' if ^ . Thus there 

should be a different stable string for each such 
fc, with string tension Uk- (Conjugate strings will 
have the same tensions.) Of course, it might be 
that we just have k fundamental strings joining 
such sources, in which case one will find ak = ka 
and we have no new string. As we shall see, this is 
not the case: we find new tightly bound fc-strings. 
The first k = 2 string appears in SU(4) and the 
first fc = 3 one in SU(6). The values of Ufe are 
interesting because they carry information about 
confinement. Also because there are specific con¬ 
jectures about their values from M (-theory) QCD 
1^, from arguments about Casimir scaling 
and from certain models, such as the bag model 
]r^ . Moreover such strings may have striking im¬ 
plications for the glueball mass spectrum as 
a function of N. 

The simplest way to calculate ak is by calcu¬ 
lating the mass of a fc-string that winds around 
a spatial torus. In a confining theory it cannot 
break and will have a length I = aL (on an 
spatial lattice). For long enough strings, the mass 
of such a loop is given by 

,,, I —^oo , 7t(D — 2) Cs /-V 

mk{l) ^ aul - - -p. (8) 

where the 0(1//) correction is the Casimir en¬ 
ergy of a periodic string and Cg is proportional 
to the central charge. This correction is univer¬ 
sal 0 since it depends only upon the massless 
modes in the effective string theory and does not 
depend upon the detailed and complicated dy¬ 
namics of the flux tube on scales comparable to 
its width. The central charge is given by the 
number of massless bosonic and fermionic modes 
that propagate along the string. In practice it 
is usually assumed that Cg = 1, corresponding 
to the simplest possible (Nambu-Goto) bosonic 
string theory. However, these modes are not re¬ 
lated to the fundamental degrees of freedom of 
our SU(N) gauge theory in any transparent way 
and the presence of fermionic modes is certainly 
not excluded. Direct numerical evidence is hard 
to get since we are interested in small corrections 
to long, massive strings. In Fig.^ we show the 
effective value of Cg that one obtains 0 by fit¬ 
ting the masses of two strings of different lengths 
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(as indicated) to eqn(|^ in SU(2) and for a lattice 
spacing a ~ ~ 0.07/to. This provides 

some evidence for the simple bosonic string cor¬ 
rection for lengths I > 1 . 2 /to. 

Assuming, then, eqn(||) with Cs = 1 we obtain 
from our calculated loop masses the string ten¬ 
sions shown in Fig.^, from which we obtain 

1.357 ±0.029 SU(4) 

a-io a \ 1.583 ±0.074 SU(5) ‘ ^ ’ 

Clearly the k = 2 string is tightly bound, and 
the string tension is incompatible, for example, 
with the naive bag model. What we find is that 
it falls between the MQCD and Casimir scaling 
predictions and, at the two sigma level, is con¬ 
sistent with both. We remark that a very recent 
higher statistics calculation ||^ favours MQCD 
over Casimir scaling. 


6. TOPOLOGY 


Gauge fields in four (suitably compactified) Eu¬ 
clidean dimensions possess non-trivial topological 
properties characterised by an integer topological 
charge Q. Topological fluctuations break the (ap¬ 
proximate) 1/^(1) symmetry of QCD, thus lead¬ 
ing to the non-Goldstone character of the 77 '. 
One can argue that if Y = 3 is close to Y = 00 
one can relate |^| the mass of the 77 ' to the topo¬ 
logical susceptibility 




( 10 ) 


of the gauge theory without quarks. (V is the 
space-time volume.) This suggests a value Xt — 
(ISOMeE)^, which is indeed close to what one 
finds in pure gauge theories ||^. All this fits to¬ 
gether if we can show that the Y = 3 values of the 
pure gauge susceptibility and of Yto^, are close 
to their Y = 00 values. While we cannot say 
much about the latter, we can and shall address 
the former. 

We calculate Q for our lattice fields by stan¬ 
dard cooling methods which for this purpose are 
reliable [|^. These calculations are performed 
simultaneously with the glueball calculations dis¬ 
cussed earlier and the values of xj'^! 6 ^" 

trapolated to the continuum limit with an c?a 



Figure 4. The D=3±l effective string correction 
coefficient estimated from the masses of flux loops 
of different lengths (indicated by the span of the 
horizontal error bar). The solid line is what one 
expects for a simple bosonic string. For compar¬ 
ison the dashed line indicates the value for the 
Neveu-Schwartz string. 
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Figure 5. The ratio of k = 2 and k = 1 string 
tensions in our D=3±l SU(4) (•) and SU(5) (o) 
lattice calculations plotted as a function of a^a. 
Extrapolations to the continuum limit, using a 
leading O(a^) correction, are displayed. 
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Figure 6. The continuum topological susceptibil¬ 
ity in units of the string tension plotted against 
\/N‘^. a linear extrapolation to iV = oo is shown. 


correction, following eqn(||). The resulting con¬ 
tinuum values are plotted against 1/iV^ in Fig.^ 
We observe that, just as for the glueball masses, 
the iV-dependence is well described by a modest 
0(1/7V^) correction and the topological suscepti¬ 
bility, expressed in physical units, does not differ 
greatly when we go from N = 3 to N = oo. 

These conclusions are somewhat weaker here 
than for the glueballa because of the rapidly grow¬ 
ing errors as N increases. The reason for this 
turns out to be interesting and easy to see ||]. 
As N grows isolated instantons become increas¬ 
ingly unlikely. This is because of the factor in the 
instanton density 

where A is the ’t Hooft coupling that we keep 
fixed as N oo. In the real vacuum this argu¬ 
ment holds for instantons with p <C 1 ‘fermi’ and 
so these instantons are exponentially suppressed 
as N increases. Now the Monte Carlo changes Q 
by an instanton shrinking through small values 
of p down to p ^ a where it can vanish through 
the lattice. (Or the reverse process.) However 


eqn(|ll|) tells us that the probability of a very 
small instanton goes rapidly to zero as N grows. 
Thus the lattice helds rapidly become constrained 
to lie in given topological sectors and for this 
quantity the Monte Carlo rapidly ceases to be 
ergodic as N grows, and the stattistical errors on 
(Q^) grow rapidly - as observed. 

We can attempt to calculate D{p) and so see 
this effect directly. To do so we examine the topo¬ 
logical charge density on the cooled fields and as¬ 
sociate an instanton to each peak in this density 
using the semiclassical formula 

Qpeak — 5 T- (12) 

This procedure is clearly of ambiguous validity 
for broad instantons, but is unambiguous for the 
high narrow peaks that correspond to the smaller 
instantons of interest to us here. In Fig.^we dis¬ 
play the result of such a calculation, comparing 
the instanton densities in SU(2), SU(3) and SU(4) 
on lattices with (almost) equal volumes and lat¬ 
tice spacings when expressed in units of the string 
tension. We clearly see a drastic suppression of 
small instantons with increasing N. 

7 . CONCLUSIONS 

We saw above that, as far as the lightest glue- 
balls, string tensions and topological susceptibil¬ 
ity are concerned, all SU(N) theories can be de¬ 
scribed by a modest leading 0{1/N‘^) correction 
to the SU(oo) limit. In this sense we can say that 
not only is TV = 3 close to N = oo but so is TV = 2. 

The large amount of interesting physics deliv¬ 
ered by these modest (workstation) calculations, 
provides a strong motivation for going further. 
Currently we |po|] are starting much larger calcu¬ 
lations on anisotropic lattices, which will enable 
us to obtain accurate mass estimates for a much 
larger range of glueball states and very accurate 
string tension ratios. Simultaneously we are 
using overlap fermions ||^ to determine the re¬ 
lationship between topology and chiral symmetry 
breaking as a function of N] and to say something 
about instantons and topology at large N. A 
more ambitious project, which would require the 
use of one of the larger ~ 1 teraflop resources that 
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Figure 7. The number density of topological 
charges plotted as a function of the charge ra¬ 
dius, p; obtained for N = 2, 3,4 from 20'^ lattices 
with Qy/a ~ 0.16. 


are becoming available to many groups, would 
be to calculate the quenched hadron spectrum at 
large N - which is interesting since it approaches 
the correct spectrum of SU(A^) QCD as we take 
N —>■ oo at fixed non-zero quark mass. And then, 
eventually, to dynamical quarks at large N and 
the 77 ' ... 

Acknowledgements 

I am very grateful to the organizers for inviting 
me to participate in this interesting (and enjoy¬ 
able) workshop, and to my fellow participants for 
many useful comments and discussions. 


REFERENCES 


1 . 


G. ’t Hooft, Nucl. Phys. B72 (1974) 461. 
E. Witten, Nucl. Phys. B160 (1979) 57. 

S. Coleman, 1979 Erice Lectures. 


A. Manohar, 1997 Les Houches Lectures, hep- 


ph/9802419 


2. M. Teper, Phys. Rev. D59 (1999) 014512. 

3. M. Teper, Phys. Lett. B397 (1997) 223. and 
unpublished. 


4. M. Wingate and S. Ohta, |hep-lat/0006016 ; 
Nucl. Phys. Proc. Suppl. 83 (2000) 381. 

5. E. Witten, Nucl. Phys. B149 (1979) 285. 

6 . B. Lucini and M. Teper, JHEP 0106 (2001) 
050. 

7. B. Lucini and M. Teper, Phys. Lett. B501 
(2001) 128; Phys.Rev. D64 (2001) 105019. 

8 . A. Hanany, M. J. Strassler and A. Zaffaroni, 
Nucl. Phys. B513 (1998) 87. 

M. J. Strassler, Nucl. Phys. Proc. Suppl. 73 
(1999) 120^ _ 

9. M. Teper, hep-th/9812187 . 

10. C. Morningstar and M. Peardon, Phys. Rev. 
D56 (1997) 4043; Phys. Rev. D60 (1999) 
034509. 

11. G. Parisi, in High Energy Physics - 1980 (AlP 
1981). 

P. Lepage and P. Mackenzie, Phys. Rev. D48 
(1993) 2250. 

P. Lepage, Schladming Lectures, hep- 


lat/9607076 . 


12. J. Ambjorn, P. Olesen and C. Peterson, Nucl. 
Phys. B240 (1984) 189, 533; B244 (1984) 262; 
Phys. Lett. B142 (1984) 410. 

S. Deldar, Phys. Rev. D62 (2000) 034509; 
JHEP 0101 (2001) 013. 

G. Bali, Phys. Rev. D62 (2000) 114503. 

V. Shevchenko and Yu. Simonov, Phys. Rev. 
Lett. 85 (2000) 1811; |hep-ph/01041^ . 

13. K. Johnson and C. B. Thorn, Phys. Rev. D13 
(1976) 1934. 

T. H. Hansson, Phys. Lett. B166 (1986) 343. 

14. M. Liischer, K. Symanzik and P. Weisz, Nucl. 
Phys. B173 (1980) 365. 

15. J. Polchinski, String Theory^ Vol 1 and II 
(CUP, 1998). 

16. L. Del Debbio, H. Panagopoulos, P. Rossi and 
E. Vicari, Ihep-th/OlOblS^ ; |hep-th/01110^ . 

17. G. ’t Hooft, Physics Reports 142 (1986) 357. 

18. E. Witten, Nucl. Phys. B156 (1979) 269. 

G. Veneziano, Nucl. Phys. B159 (1979) 213. 

19. M. Teper, Nucl. Phys. Proc. Suppl. 83 (2000) 
146. 


20. B. Lucini, M. Teper and U. Wenger, in 
progress. 

21. N. Gundy, M. Teper and U. Wenger, in 
progress. 

22. H. Neuberger, these proceedings. 
























